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The critical current of the spin transfer torque-driven magnetization dynamics was studied by 
taking into account both spin pumping and the finite penetration depth of the transverse spin 
current. We successfully reproduced the recent experimental results obtained by Chen et al. [Phys. 
Rev. B 74, 144408 (2006)] and found that the critical current remains finite even in the zero- 
thickness limit of the free layer. We showed that the remaining value of the critical current is 
determined mainly by spin pumping. We also showed that we could control the critical current by 
varying the spin diffusion length of the nonmagnetic electrode adjacent to the free layer. 
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Spin transfer torque (STT) is the torque due to the 
transfer of transverse spin angular momentum from the 
conducting electrons to the magnetization of a ferromag- 
net P, 0|. The STT in magnetic multilayers such as 
the current perpendicular-to-the-plane giant magneto- 
resistive (CPP-GMR) a i, n and tunnel magneto- 
reisistive (TMR) 0, 0, [1| spin valves has been investi- 
gated intensively because STT-driven magnetization dy- 
namics is a promising technique to operate the spin- 
electronics devices such as magnetic random access mem- 
ories and microwave oscillators. One of the main ob- 
stacles in developing STT-based spin-electronics devices 
is the high critical current density. The critical current 
density required to induce the STT driven magnetization 
dynamics in CPP-GMR spin valves is as high as 10 6 — 10 8 

A/cm 2 0, Ea, E3, E3 - 

On the other hand, the CPP-GMR spin-valve is one 
of the promising candidates for the read head for ultra- 
high-density magnetic recording pH , [Til ], ft is known 
that STT-driven magnetization dynamics produces noise, 
and that low critical current density is required for the 
read head application [15|. Therefore, it is natural to ask 
how to control the critical current density of STT-driven 
magnetization dynamics in CPP-GMR spin valves. 

STT was first proposed by Slonczewski [l[ and inde- 
pendently by Berger in 1996. In Slonczewski's theory 
the critical current of STT-driven magnetization dynam- 
ics is expressed as pjl Qj} , 



2eMSd 
I c = — z a uj, 



(1) 



where e is the elementary charge and h is the Dirac con- 
stant. M, S, d, 7 and ag are the magnetization, cross- 
section area, thickness, gyromagnetic ratio and intrinsic 
Gilbert damping constant of the free layer, respectively. 
w is the angular frequency of the magnetization around 
the equilibrium point. The transverse spin polarization 
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FIG. 1: The nonmagnetic (N) / ferromagnetic(F) multilayer, 
i.e., CPP-GMR spin valve we consider is schematically shown. 
The symbols are defined in the text. 



coefficient r\ depends only on the relative an gle of the 
magnetizations of the fixed and free layers [TL 1 1 6j] . Ac- 
cording to Slonczewski's theory, we can control the crit- 
ical current by varying the thickness of the free layer, d 
and the critical current vanishes in the limit of d — > 0. 

However, recently, Chen et al. [l2j reported that the 
critical current of STT-driven magnetization dynamics 
of a CPP-GMR spin valve remains finite even in the 
zero-thickness limit of the free layer. What are missing 
from the above consideration based on Slonczewski's the- 
ory are the effects of the finite penetration depth of the 
transverse s pin current, At, [13, llSl l2fl| an d spin pumping 
HH HI IM IM H - The Penetration depth of the trans- 
verse spin current is the characteristic length of the fer- 
romagnetic metal over which the transfer of the spin an- 
gular momentum from conducting electrons to the mag- 
netization is achieved. If the free layer is thinner than 
A t , the conducting electrons cannot transfer their angu- 
lar momentum to the magnetization to exert STT. Spin 
pumping is the phenomenon by which the spin current is 
pumped out from the free layer into the other layers. The 
magnetic (Gilbert) damping of the free layer is enhanced 
by spin pumping. Therefore, we need to analyze the ex- 
perimental results by taking into account both the finite 
penetration depth of the transverse spin current and spin 
pumping to understand the mechanism that determines 
the critical current of STT-driven magnetization dynam- 
ics in magnetic multilayers. 

In this paper, we study the critical current of STT- 
driven magnetization dynamics by taking into account 
both the finite penetration depth of the transverse spin 
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current and spin pumping. In order to analyze the exper- 
iments by Chen et al. [12|, we extend the spin pumping 
theory with the finite penetration depth [201 ] to include 
the electric current. We show that the critical current 
remains finite even in the zero-thickness limit of the free 
layer, which agrees quantitatively well with the results 
of Ref . fl2| , and that the remaining value is determined 
mainly by spin pumping. We find that we can control the 
remaining value of the critical current by varying the spin 
diffusion length of the nonmagnetic electrode adjacent to 
the free layer. The longer the spin diffusion length of the 
nonmagnetic electrode, the smaller the remaining value 
of the critical current. 

The system we consider is shown in Fig. [TJ Two fer- 
romagnetic layers (Fi and F2) are sandwiched by the 
nonmagnetic layers Nj (i = l — 7). The Fi and F2 layers 
correspond to the free and fixed layer, respectively, 
(k = 1, 2) is the unit vector pointing the direction of the 
magnetization of the F& layer, d k and Lj are the thick- 
nesses of the Ffe and Nj layers, respectively. The electric 
current I flows from the N7 layer to the Ni layer. 

In order to analyze the STT-driven magnetization dy- 
namics in the multilayer system shown in Fig. [TJ we ex- 
tend t he sp in-pumping theory with the finite penetration 
depth [20] to include the electric current. The electric 
current and pumped spin current at the Ffc/Nj interface 
(into Nj) are derived by the circuit theory (25[, and ex- 
pressed in terms of the charge accumulation /XN;,F fc and 
the spin accumulation fi^ i p k as [23l. [25j 

/Ffc/Nl = l {2(AiFfc -^Nj+pm fe .( MFfc - MNi )}, (2) 
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where h = 2irh is the Planck constant, g = g' ' + g 
is the sum of the spin- up and spin-down conductances, 
P = {9 9 )/(g +g ) is the spin polarization of the 
conductances and g T n-\ is the real (imaginary) part of the 
mixing conductance. The spin current at each Ffc/Nj and 
Nj/Nj interface (into Nj) are given by [19, [2^] 
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where tjix is the real (imaginary) part of the transmission 

mixing conductance at the Ffc/Nj interface and gNi/Nj is 
the conductance of the one spin channel at the interface. 
The spin current of Eq. §4§ is obtained from the circuit 
theory of Brataas et al. [251 ] eliminating the assumption 
that the non-equilibrium distribution function of the elec- 
trons in a ferromagnetic layer is aligned to the direction 



of the magnetization in spin space. It should be noted 
that the transmission mixing conductance in Eq. fll]) is 
different from that defined by Zwierzycki et al. [261 ]. 
Zwierzycki et al. calculated a transmission mixing con- 
ductance defined through a N/F/N junction defined as 
t 1l = t r t l* where t a = t 
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the tranmission coefficient for electrons from F (N) to N 
(F), and showed that depends on the thickness of the 
ferromagnetic layer d due to the phase factor e ikj - d [26| . 
On the other hand, the transmission mixing conductance 
in Eq. (g]) is defined by iLL = Re(Im)[t] r ^ N tp*_^ N ], and 
independent of the thickness of the ferromagnetic layer. 
Although the original formulation of the circuit theory 
assumed the spatially uniform charge and spin accumu- 
lation [251] . it has been shown that the circuit theory is 
applicable to the diffusive system [24], [53] ■ It should be 
noted that there is a controversial issue regarding the 
transverse spin accumulation in the ferro mag netic layer, 

i4 = m x (/x F x m) [H M, M, HE HI m El. 

The spin accumulation in the nonmagnetic layer, /lin 
obeys the diffusion equation [3l|, and is expressed as a 
linear combination of exp(±a;/A s d(N)), where A s( j(n) is the 
spin diffusion length of the nonmagnetic layer. The spin 
current in the nonmagnetic layer is given by 



dx 2e 2 



(6) 



where on is the conductivity of the nonmagnetic layer. 

The longitudinal spin accumulation in the ferromag- 
netic layer, /Xp = (m fiF)m, also satisfies the diffusion 
equation, and is expressed as a linear combination of 
exp(±x/A s d(F L ))j where A s( j(f l ) is the lon gitu dinal spin 
diffusion length of the ferromagnetic layer [3l| . The lon- 
gitudinal spin current in the ferromagnetic layer is 

(m-rj)m= -A ^(^J, -<j F n F )m, (7) 

where fj,p ^=J E ^def^^ andcrp ^ are the electro-chemical 
potential and the conductivity for the spin-up (spin- 
down) electrons, respectively. The spin polarization of 
the conductivity is defined as = (cr| — erp ) / (op + erp ) . 

The transverse spin accumulation in the ferromagnetic 
layer obeys [HJ 



dx 2 
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where Xj = y (D F + D F )H/ (2 J) and A sd (p T ) is the trans- 
verse spin diffusion length. J is the strength of the ex- 
change field [l9j] and D F ^ is the diffusion constant of 
spin-up (spin-down) electrons. The spin polarization of 
the diffusion constant is defined as /3'=(D F — Dp) / (D F + 
Dp). The transverse spin accumulation is expressed as 
a linear combination of exp(±cc//+) and exp(±a;/Z_), 



where l/l± = J (l/A 2 d( . FT - ) )=F(i/Aj). The penetration 
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depth of the transverse spin current At is defined as 
1/At = Re[l/Z + ] [20j|. The transverse spin current in the 
ferromagnetic layer is expressed as 
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where 4^ = [4/(1+^+4/(1-^)1/2 [H,[I 



The total spin currents across the N3/F1 and F1/N4 in- 
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terfaces, i.e., I, 

exert the torque r = mi x [(Is ±; + Is ZJ ) x mi] on the mag- 

f 1 2) 

netization mi . In order to obtain the spin current I« , 
we solve the diffusion equations of spin accumulations in 
each layer. The boundary conditions are as follows. We 
assume that the thickness of the Ni and N7 layer, L\ 
and -L7, are sufficiently thick enough compared to their 
spin diffusion length, and that the spin current is zero at 
the outer boundary of the Ni and N7 layer. We also as- 
sume that the spin current is continuous at all F/j /Nj and 
Nj /~Nj interfaces and that the electric current is constant 

(1 2) 

through the entire structure. The spin current I, is 
obtained as a function of the electric current I and the 
pumped spin current IP ump . 

The torque r modifies the Landau-Lifshitz-Gilbcrt 
(LLG) equation of the magnetization mi. The LLG 
equation conserves the magnitude of the magnetization, 
and thus the vectors rhi and mixrhi are perpendicular 
to the magnetization mi. Since the torque r is perpen- 
dicular to mi the torque can be decomposed into the 
directions of rhi and rnjXriii. The LLG equation of mi 
is expressed as [2(| [H, [HJ 
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(10) 



where B c ff is the effective magnetic field, a' = a c +a pump 
represents the enhancement of the Gilbert damping con- 
stant. The enhancement a c is proportional to the electric 
current / and independent of the pumped spin current 
jpump rj-^g gnhancgujgnf; a pump represents the contribu- 
tion from the pumped spin current and is independent of 
the electric current. It should be noted that the enhance- 
ment a pum p differs from the result of the conventional 



spin-pumping theory [24J because a pump is a function of 



At. The enhancement of the gyromagnetic ratio 7 e ff/7 
is also a function of the electric current and the pumped 
spin current. 

Let us move to the analysis of experimental results of 
Ref. [12}. In general, the dynamics of the magnetiza- 
tion mi determined by Eq. (jTUJ) is very complicated; 
thus, wc cannot obtain the analytical expression of the 
critical current of STT-driven magnetization dynamics 
of the magnetization mi. However, in the experiment 
of Ref. [l2| , the system, and therefore the dynamics of 
mi, have axial symmetry along the direction normal to 
the film plane because the high magnetic field (about 7 



T) is applied along this direction. Then we assume that 
the magnetization of the Fi layer mi precesses around 
the magnetization of the F2 layer rri2 with the relative 
angle 9 and the angular frequency u). The critical current 
of STT-driven magnetization dynamics is defined by the 
current that satisfies the condition, olq + a c + cx puiap = 0. 
The critical current I c is expressed as 



2eM3d 1 . . 
= — ftrfij — + a pumpM 



(11) 



where r) is the effective transverse spin polarization co- 
efficient that is determined by the diffusion equations of 
the spin accumulations, and thus the coefficient fj is the 
function of di/A S( j(F L ) an d d 1 /l±. 

The parameters we used are as follows. The system 
consists of nine layers shown in Fig. [TJ where Fi and 
F 2 are Co, Ni, N 3 , N 4 , N 5 and N 7 are Cu, and N 2 and 
Ng are Pt. The thicknesses of the N3, N4 and N5 layers 
are 10 nm, the thicknesses of the N 2 and N$ layers are 
3 nm and the thickness of the F 2 layer is 12 nm [l2l ]. 
The thickness of the Ni and N7 layers are taken to be 
10 /zm, which is sufficiently longer than the spin diffu- 
sion length. The resistivity (^n)" 1 of Cu and Pt are 14 
and 42 fmm, respectively [33|. The spin diffusion length 
An of Cu and Pt are 1000 and 14 nm, respectively [331 ] . 
The conductance at the Cu/Pt interface gcu/vt/S is 35 
nm~ 2 [33]. The magnetization M, the intrinsic Gilbert 
damping constant ao and the gyromagnetic ratio 7 of 
Co are 0.14 T, 0.008 and 1.89 x 10 11 Hz/T, respectively 
[HI, [13] ■ For simplicity, we assume that p=(3=(3 =0.46 for 

Co [HI]. The resistivity (4+4) _1 ano - tne longitudinal 
spin diffusion length A sd (p L ) of the Co are 60 Slnm and 
40 nm, respectively [331 ]. The transverse spin diffusion 
length is A ao v FT ) = A sd ( FL )/y / l — fi 2 [HI]. Xj is taken to 
be 3.0 nm [19[, i.e., At =4.2 nm. The conductances at the 
Co/Cu interface, g/S, qj l IS and q ; n /S, are 50, 27 and 
0.4 nm- 2 , respectively ||I3, EElSE EH]- The angular 
frequency is uj — 7(-B apP i + 47rM) where the strength of 
the applied magnetic field £? app i is 7 T [ijj] ■ The relative 
angle of the two magnetizations 9 is assumed to be 0.997T 
[12j . Although there are many material parameters in 
our calculation these values except t t \ are determined by 
the experiments and ab initio calculations. The value of 
tlyS is determied by fitting, and taken to be 6.0 nm~ 2 . 
According to Ref. 12], the experimental values are the 
low temperature values. 

The obtained critical current density is plotted by a 
solid line against the thickness of the free layer d\ in Fig. 
[2] The experimental results of Ref. [l2j] are shown by 
open circles. One can see that our results agree well with 
the experimental results. The critical current density de- 
creases as the thickness of the free layer decreases, and 
remains finite even in the zero-thickness limit of the free 
layer. In order to see the main mechanism that deter- 
mines the remaining value of the critical current density, 
we decompose J c of Eq. (fTTj) into two parts as I c = 1®+!%, 
where 7° is the component proportional to ao and I£ is 
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FIG. 2: The critical current density vs the thickness of the 
free (Fi) layer. The circles are the experimental result of 
Chen et al. for 70x140 nm 2 junctions The solid line 

corresponds to I c /S (see Eq. JHJ-) The dotted line and 
dashed-dotted line correspond to I® / S and Ic/S, respectively. 
The dashed line corresponds to I c /S in the limit of A t — > 0. 

the component proportional to a pump . In Fig. [2l the 
components 1° / S and I£/S are plotted by dotted and 
dot-dashed lines, respectively. As shown in Fig. [5]the re- 
maining value of the critical current in the limit of d\ — ► 
is determined mainly by the spin pumping. Although 
I®/S is also finite in the limit of d\ — > because of the 
finite penetration depth of the transverse spin current At 
in the Fi layer, the remaining value is small compared to 
Ic/S. The dashed line in Fig. [2] is the calculated critical 
current I c /S in the limit of A t — > 0. According to Fig. [3] 
we conclude that the effect of the finite penetration depth 
At is less important to describe the results of Ref . [12| . 

The reason why both 1^ and If remain finite in the 
limit of di — » is understood as follows. Slonczewski as- 
sumed that the transverse spin current injected into the 
free layer is absorbed at the interface, and thus, STT is 
independent of the thickness of the free layer. The critical 
current is determined by the competition between STT 
and the magnetic (Gilbert) damping of the free layer. 
The spin relaxation due to the Gilbert damping is pro- 
portional to the thickness of the free layer di, and thus, 
the critical current given by Eq. (Q]) is proportional to 
d\ and vanishes in the limit of d\ — > 0. If the penetra- 
tion depth of the transverse spin current At is finite, the 
transverse spin current is not fully absorbed in the free 
layer in the case of d\ <C At . Then the strength of STT is 
decreased compared to the prediction of Slonczewski [l[ , 
and thus, the critical current is increased. Spin pumping 
enhances the Gilbert damping, and the spin relaxation 
due to spin pumping is independent of the thickness of 
the free layer. Thus, JP remains finite in the limit of 
di -> 0. 

The dependences of the remaining value, about 1.6 x 
10 8 A/cm 2 , on the parameters given above are as follows. 
If the resistivity and the longitudinal spin diffusion length 
of Co are taken to be 210 fmm and 38 nm, respectively, 
which are the room temperature values [33[ , the remain- 
ing value is estimated to be 1.5 x 10 8 A/cm 2 . The reduc- 
tion the longitudinal spin diffusion length decreases the 
penetration depth At, and thus, the remaining value is re- 



duced. The values of conductances, g and g±\, include the 
effect of the Sharvin conductace [26j |. If g/S, gj^/S and 
<7j / 'S are taken to be 19.3, 14.6 and -1.1 nm -2 , respec- 
tively, which are the bare values estimated by ab initio 
calculation (26[, the remaining value is estimated to be 
1.0 x 10 8 A/cm 2 . The reduction of the mixing conduc- 
tance decreases the effect of spin pumping [231 ] . and thus, 
the remaining value is reduced. If the transmission mix- 
ing conductance t].\/S is taken to be 3.0 (12.0) nm -2 , 
which is half (twice) compared to the value used in Fig. 
[SJ the remaining value is estimated to be 1.5 (1.8) x 10 8 
A/cm 2 . The reduction (enhancement) of the transmis- 
sion mixing conductance decreases (increases) the effect 
of the transverse spin accumulation, or equivallently the 
penetration depth, on the spin current given by Eq. 
and thus, the remaining value is reduced (enhanced). We 
conclude that although there are many material param- 
eters in our calculation the parameter dependence of the 
remaining value is small, and our calculation gives the 
correct order of the critical current. 

The above results imply that we can increase or de- 
crease the critical current by controlling the spin pump- 
ing. Spin pumping is the phenomenon by which the pre- 
cessing magnetization of the free layer pumps spin cur- 
rent into the other layers. The other layers act as an 
additional spin sink and the magnetic damping of the 
free layer is enhanced by spin pumping. The ability of 
the spin sink is determined by the spin diffusion length 
since the spin diffusion length is inversely proportional 
to the square root of the spin scattering rate. Materi- 
als with short (long) spin diffusion length act as a good 
(bad) spin sink. One may expect that if the nonmagnetic 
layer adjacent to the free layer is made of material with a 
long spin diffusion length, the Gilbert damping constant 
and therefore the critical current is suppressed. In the 
limit of infinite spin diffusion length, An — ► oo, there is 
no spin flip scattering in the nonmagnetic layer and the 
spin pumping into the nonmagnetic layer is forbidden. 
On the other hand, if the nonmagnetic layer adjacent to 
the free layer is made of a material with a short spin 
diffusion length, the Gilbert damping constant and the 
critical current is enhanced. In the limit of An —> 0, the 
pumped spin current is absorbed at the interface and en- 
hancement of the critical current due to spin current is 
maximized. 

In order to verify the above statement, we an- 
alyzed the critical current of the five-layer system, 
N1/F1/N4/F2/N7 (see Fig. Q]), where all parameters ex- 
cept the spin diffusion length of the Ni layer, Ani , are the 
same as those used in the analysis of Chen's experiment. 
In Fig. [3l we plot the critical current in the zero-thickness 
limit of the free layer as a function of the spin diffusion 
length of the N4 layer An x ■ One can see that the critical 
current is a decreasing function of An x ■ The critical cur- 
rent remains finite in the limit of An — > 00 because of the 
spin pumping into the N4 layer and the finite penetration 
depth of the transverse spin current. The result shown 
in Fig. [3] shows that we can control the critical current 
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FIG. 3: The dependence of the critical current in the zero- 
thickness limit of the free layer in the N1/F1/N4/F2/N7 five- 
layer system on the spin diffusion length of the Ni layer, Anj ■ 

by varying the spin diffusion length of the nonmagnetic 
electrode adjacent to the free layer. 

We cannot apply the present formula directly to a 
magnetic tunnel junction (MTJ) because spin accumu- 
lation is not well-defined in an insulator (I). Although 
the spin pumping across the insulating barrier is beyond 
the scope of this paper, the spin pumping into the metal- 
lic electrode should give the finite remaining value of 
the critical current. Recently spin pumping in a mag- 
netic tunnel junction (MTJ) was studied by Moriyama 
et al. [36]. They studied a ferromagnetic resonance in 



Al/A10/Ni 80 Fe 2 o/Cu MTJ, and found the generation of 
the voltage on the order of few /xV, which is qualitatively 
explained by the theory of spin pumping in a metallic 



structure |37j, but requires unreasonably large value of 
mixing conductance. The results of Moriyama et al. 
suggest that a new non-equilibtirum phenomena exists 
in MTJ, e.g. charge pumping or the development of the 
longitudinal spin accumulation in a ferromagnetic layer 

In conclusion, we studied the critical current of spin 
transfer torque-driven magnetization dynamics by tak- 
ing into account both the finite penetration depth of the 
transverse spin current in the ferromagnetic layer and 
spin pumping. We extend the spin-pumping theory with 
the finite penetration depth to include the electric cur- 
rent and successfully reproduced the experimental results 
of Ref . [12| . We showed that the critical current remains 
finite in the zero-thickness limit of the free layer and the 
remaining value is determined mainly by spin pumping. 
We also showed we can control the remaining value of 
the critical current by varying the spin diffusion length 
of the nonmagnetic electrode adjacent to the free layer. 
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